Hubbard-I approach to the Mott transition 
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A Hubbard-I like approach is reconsidered in the context of the metal-insulator transition in 
strongly correlated electron systems. We show that this approach yields a simple treatment of the 
transition (at half-filling) and the metallic phase (for arbitrary fillings) in a constrained Hubbard 
model, that is equivalent to a particular case of a Hubbard bond charge Hamiltonian on bipartite 
lattices. These calculations, which are to our knowledge the first Hubbard-I calculations for this 
model, yield results which compare well with known results for this model. For half filling, we also 
compare these results with those obtained using the Gutzwiller approximation. Finally, we propose 
a candidate variational wave function connecting the Gutzwiller wave function approach with the 
scheme presented here. 
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Introduction The metal to Mott insulator (MI) tran- 
sition, envisaged by Mott [1], is one of the striking ef- 
fects induced by strong electronic correlations in many 
electron systems. The Coulomb interaction between elec- 
trons U leads to highly correlated ground states in these 
systems, rendering their description quite challenging. 
As a result there have been various theoretical attempts 
at providing a satisfying description of the Mott transi- 
tion. 

The first of these was due to Hubbard [2], who pro- 
vided a seminal nonperturbative approach - the so-called 
Hubbard-I (H-I) approximation - to a simplified interact- 
ing electron problem described by the Hubbard model. 
H-I describes both (i) the atomic limit (i.e. limit of van- 
ishing bandwidth W — > 0) exactly, in particular yielding 
two atomic levels corresponding to single and double lo- 
cal occupancy, and (ii) the non- interacting case (U = 0) 
exactly, and so held some promise of describing the in- 
termediate physics in a consistent interpolating fashion. 
For finite bandwidth, the atomic levels broaden into two 
"dynamic" (sub)bands which are occupation-number and 
interaction dependent. These are always split by a gap 
for all C/ > indicating an insulator phase. Unfortu- 
nately, a metal-insulator transition is predicted exactly 
at U = 0, associated with gap closure. The situation was 
somewhat improved in the Hubbard-III approximation 
[3], where scattering and resonance broadening correc- 
tions shift the transition point to U c ~ W. However, 
this approach does not predict the expected Fermi Liq- 
uid properties on the metallic side, [4] . These drawbacks 
not withstanding, the Hubbard approaches are very im- 
portant for the conceptual introduction of the Hubbard 
subbands and also because they give a basic description 
of the insulating phase. 

A complementary approach, starting from the weak 
coupling limit, is based on the Gutzwiller variational 
wave function [5], which describes an increasingly cor- 
related Fermi liquid with increasing interaction U. 
Brinkman and Rice [6] showed, using the so-called 



Gutzwiller approximation [5], that increase in U is asso- 
ciated with the diminishing of the quasiparticle residue 
of the Fermi liquid, or equivalcntly with the increase 
in quasiparticle effective mass. In this framework, a 
metal-insulator transition occurs when the effective mass 
diverges and is thus driven by quasiparticle localisa- 
tion. The critical value U c obtained in the (Gutzwiller)- 
Brinkman-Rice (GBR) approach is much larger than in 
the Hubbard approach. This essentially weak coupling 
method gives a good low energy description of the metal, 
but does not describe the precursors of the Hubbard 
bands which should plausibly appear on the metallic side. 

More involved contemporary methods, such as dynam- 
ical mean field theory (DMFT) [7] or Green's function 
ladder diagram techniques [8], provide a connection be- 
tween a strongly correlated fermi liquid picture on the 
one hand, and formation of Hubbard bands on the other. 
While providing greater insight into the physics at play, 
these methods do not provide a simple picture character- 
istic of more analytical techniques. In this paper, we re- 
consider the Hubbard-I approach as a systematic means 
of obtaining a picture, complementary to the one ob- 
tained using the Gutzwiller approximation, of the Mott 
transition from the point of view of strong coupling the- 
ory. We shall show that there are strongly correlated 
Hamiltonians for which the Hubbard-I approach yields 
good qualitative and quantitative results. 

Inspection of the Hubbard-I approach Contrary to 
early expectations, the Hubbard approach, based on a 
Green's function decoupling scheme, is now seen as a 
large-U approximation, despite its non perturbative na- 
ture. One may therefore suspect that some problems 
encountered in the Hubbard approach may originate in 
limitations of strong-coupling perturbation theory. 

Consider therefore the Hubbard model: H = T + 
UD, where D = n^n^ is the number of dou- 
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bly occupied sites and t = J2i j a 
netic energy. The large-U perturbational expansion 
starts from splitting the kinetic energy into three parts 



T = T + T+i + T_i. T = T UH b + T LHB is the 
double occupancy conserving hopping (commuting per- 
turbation) in the Upper and Lower Hubbard Bands 

nia-)c ia .Cj a (l — nja-), corresponding to hopping of pro- 
jected electrons on doubly and singly occupied sites re- 
spectively, and interband hopping terms T+i and T_i = 
TL (off-diagonal perturbation). The perturbational ex- 
pansion for the Hubbard model is well known and can be 
performed e.g. using the method of canonical transfor- 
mations [9, 10] which to second order yields the effective 
Hamiltonian: 

H ( §=% + ub+±j[f +ll f_ 1 ). (1) 

Recall that the perturbation expansion to any or- 
der leads to an effective Hamiltonian eliminating mix- 
ing between the degenerate subspaces of the unperturbed 
Hamiltonian U. Due to this property, the average of op- 
erators jointly changing the total number of doubly oc- 
cupied sites, such as e.g. {cj a nj S c\ a hi^) ,(hj& = 1 — n,j S 
is the on-site hole occupation number), arc identically 
zero leading to the vanishing of the associated Green's 
functions. Thus in this framework, the single particle 
Green's function G(j,l) = {(cj a \c' la )) decomposes into 
a sum of two Green's functions Tf^ = ((cj a nja\ci a .nia}) 
and rj ; = {{cj a hj S \c\ a hi s )) related to the propagation of 
fcrmionic quasiparticles, in the upper and lower Hubbard 
bands. 

Consider the Hubbard-I approach to the simplest non- 
trivial level of strong-coupling perturbational expansion, 
described by the first two terms in the expansion Eq.(l) 

H c = f + UD. (2) 

In the frequency domain, the equations of 
motion for the basic Green's function are 
LoT^iu) = ^(ntg-)Sij + (([Q .n i5 .,iJ c ]|c] cr n i5 .)) u and 

<Jrfj(u) = ±(h l9 )5 }j + (([ctoh^H^hje))*,. We 
perform the following type of decoupling on the 
higher order Green's function in the equations for 

({hl+5aCl +s ^ClaCla\& jo .)) u « {c is Ci a ) ((cj' +fo /l Z+5(T |c] (T )) w , 

with analogous complementary treatment of the set 
Tfj(u)). The last two decouplings lead to magnetic and 
superconducting order parameters are set here to zero, 
as we are interested in the description of a pure Mott 
transition. Solving the equations of motion, one obtains 
the following momentum (or Fourier transformed) 
Green's functions for the system: 



The electron Green's function thus has the character- 
istic Hubbard two pole form describing two subbands, 
however the poles here have an elementary form in com- 
parison with the solution for the full Hubbard model [2] . 
The main result stemming from these Green's functions 
is that a Mott transition is elevated to finite interaction 
U. Indeed, for a paramagnetic configuration at half fill- 
ing (n-f) = (n-j-) =1/2 the lowest energy in the UHB is 
U - W/4 while the highest energy in the LHB is W/4 
(where W is the free electron bandwidth), hence a gap 
opens at a critical value U c = W/2. A fundamental 
drawback of the Hubbard-I approach to the full Hub- 
bard model is that the resultant Green's functions do not 
conserve the electron-hole symmetry at half-filling. Im- 
portantly, the presented solutions Eqs.(3,4) are explicitly 
electron hole symmetric and do not suffer from this draw- 
back. Accordingly, we show below that the gap opening 
at U = W/2 is associated with a Mott transition only 
exactly at half-filling. These results provide explicit evi- 
dence that the source of problems in the Hubbard-I ap- 
proach stem from the effects of the T±i terms. 

The Mott transition The constrained Hamiltonian H c 
Eq.(2) is interesting in itself because it can be viewed 
as a particular model of extremely correlated electrons. 
Recall that the term extremely correlated electron sys- 
tems has recently been introduced by Shastry [11], em- 
phasizing the appearance of noncanonical fcrmions re- 
sulting from the prohibition of double occupancies in the 
U = oo limit, which is a special case of symmetry of 
conserved double occupancies. The model H c considered 
here, that allows for double occupancies which are con- 
served and thus describes hopping of noncanonical (con- 
strained) fcrmions Cj a rijg and Cj a hj S can therefore be 
considered a generalization of the problem of extremely 
correlated fcrmions for general interaction strength U . 
Note that the model H c contains the U — > oo limit of the 
Hubbard model. Even for small U and low fillings, we 
discuss below that the ground states of H c correspond to 
the ground states of the U = oo Hubbard model. 

We now consider the details of the Hubbard-I descrip- 
tion of the zero temperature Mott transition in the model 
H c of extremely correlated electrons. The analysis is car- 
ried out only for the paramagnetic phases, with = nj_. 

From the band structure of the one particle Green's 
function, it is evident that the insulator phase exists only 
if U > U c and the lower band is completely filled while 
the upper band is empty. The Green's functions Eqs.(3,4) 
lead to the following equation for the number of electrons 
of a given spin species: 

no- = (1 - n s ) I f((l - n s )e - fi)p(e) de 
J-W 

+n a / f(n ff e + U - fx)p(e) de (5) 
where /(. . .) is the Fermi-Dirac distribution and p(e) is 
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the density of states. The boundaries of the insulator 
phase are obtained, independently of the lattice, when 
either \i = (1 — n s )W/2 is at the end of the lower band or 
fi = —n s W/2 + U is at the beginning of the upper band. 
For both these cases, at zero temperature, Eq.(5) reduces 
to Tifj = 1 — rig, showing that the transition only occurs at 
half-filling. The jump in p at half-filling n-j- = nj. = 1/2 
reduces to zero at U = W/2 indicating the transition 
point. The Mott "lobe" is shown in Fig.l. Note that the 
kinetic energy and interaction energy are both zero in the 
Mott phase, which is an exact feature of the model. 

Outside the Mott phase, the system is in a metallic 
phase described by four (two per spin species) dispersion- 
less extremely correlated Fermi liquids (ECFLs) associ- 
ated with the lower and upper Hubbard bands. Depend- 
ing on the filling factor and value of U, for each spin 
species one can obtain a situation where only one kind 
of Fermi liquid appears to which there pertains a single 
well defined Fermi surface (see Fig. 1). For n < 1, this 
case physically resembles the ground state properties of 
the infinite U limit of the Hubbard model. The single 
Fermi liquid comes about because for small n.f,n± the 
lower Hubbard band is wide while the upper Hubbard 
band is very narrow. Even for U < W/2, although the 
two bands overlap, the upper band centered around U 
is too narrow to be occupied in the ground state. On 
the other hand, closer to half filling n < 1 for U < W/2 
when the upper band can also be occupied, two coexis- 
tant Fermi liquids emerge, per spin species, associated 
with two Fermi surfaces. The boundary between these 
two cases is a Fermi surface topology changing transition 
of the Lifshitz type. 

We now focus on the interaction driven Mott transition 
from the metallic phase. Unlike in the Gutwziller or slave 
boson methods, the hopping in the approach discussed 
here cannot be used as an indicator of the transition, as 
it remains constant. The density of doubly occupied sites 
D = (D)/N (occupation of the upper Hubbard band) is 
a good parameter at half filling capturing the transition 
for this model. It takes a particularly simple form 

D={D)/N= \ I U p(e)de (6) 

Z J -W/2 

and depends only on U and the density of states. 
We show results for a few different density of states 
(DOSes)in Fig.2. 

Note that an important feature of our Hubbard-I cal- 
culation is that the critical point does not depend on 
the lattice type. Furthermore, lattice dependent proper- 
ties emerge for the number of doubly occupied sites. A 
linear dependence is associated here only with a rectan- 
gular density of states. In fact, the critical behaviour 
of D has universal properties depending only on spa- 
tial dimensions of the lattice. It is governed by the 
behaviour of the DOS at the bottom of the band. In- 
deed for parabolic energy profiles near the bottom of 
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Figure 1: Ground state phase diagram in the U — n plane. 
In all calculations we assume only nearest neighbour hopping 
tij — —t. A metallic ECFL exists for all n < 1. At low filling, 
the metal resembles the U = oo with no double occupancies. 
Closer to n = 1, the metal has double occupancies and is 
described by two Fermi liquids per spin species. In the H- 
I calculation, the boundary between these two behaviours is 
a Lifshitz-type transition, shown for hypercubic lattices in 
1,2,3 dimensions and for the rectangular DOS. Shown is also 
the exact result for Id, from [17] with rescaled U — > U/2, 
for comparison. All transition lines cross the n = 1 line at 
finite U — W/2 marking the point of Mott transition. Solid 
line depicts the Mott phase. The n > 1 can be obtained by 
invoking electron- hole symmetry. Inset: Ground state phase 
diagram in the U — p plane. The electron hole symmetry line 
(dotted) separates n < 1 and n > 1 cases. The n = 1 case 
coincides with this line for U < W/2 whereas for U > W/2 
opens into a Mott insulator "lobe" of finite width. 



the band, as £ 4 — W/2 => Ae^ oc fc 2 , we have 
for 1 dimension p{— W/2 + Ae) oc 1/yAe, 2 dimen- 
sions that p(— W/2 + Ae) oc const., and 3 dimensions 
p(-W/2 + Ae) oc s/Ae. Then we obtain that in Id: 
D oc [(U c - U)/U c ]i , 2d: Dot (U c - U)/U c while in 3d: 
Doc[(U c -U)/U c ]%. 

Relation to the Bond-charge Hubbard model While 
the results presented in the previous paragraph provide 
a clear and simple quantitative picture of a paramagnetic 
transition as well as the metallic phase, it is not a pri- 
ori clear how good an approximation the Hubbard-I ap- 
proach yields for the constrained model H c . This model is 
related to the class of generalized Hubbard models, dif- 
fering from the Hubbard Hamiltonian by an additional 
bond charge interaction term (see e.g. [13-15]): 



-^bond — charge T -\- U D 



■E 



Xij {jlia 



' l 3cr) c icr C iv 



(7) 



At a symmetry point, when the bond charge interaction 
is equal to the hopping = —tij, this generalized Hub- 
bard model conserves the number of double occupancies. 
In fact, the resultant symmetric model can be mapped 
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Figure 2: The double occupancy per site D at half- filling 
n = 1 as a function of U. Shown are results for hypercubic 
lattices in 1, 2, 3 dimensions and for rectangular DOS. For 
all DOSes D = 1/4, at U = 0, as in the uncorrelated metal, 
although it is a ECFL here. The Mott transition takes place 
at U = W/2. Note different critical behavior in 1,2 and 3 d. 
The curve calculated for Id coincides with known exact result 
[17] for rescaled U -> (7/2. 

onto the model H c on bipartite lattices, via the canonical 
transform U = exp[— mtJ^- Rj n jt n jl\> where Rj = ±1 
on different sublattices, as shown in [16]. Analytical 
ground states of the symmetric bond-charge model were 
obtained in [15] in the regime U > 2W at half filling, 
for arbitrary dimensions d. These results were improved 
upon in [17] revealing that a metal-insulator transition 
occurs at U = W, probably with ferromagnetic polariza- 
tion on the metallic side (for d > 1), which is supported 
by numerical evidence [18]. 

The special case of d = 1, was exactly solved in [17] 
leading to the following main results: the metal to in- 
sulator transition occurs at U c = W, and also the num- 
ber of double occupancies is D = l/27rArcCos([//[/ c ). 
In d = 1, the H-I approach used here underestimates 
the transition point, but remarkably the average num- 
ber of double occupancies calculated for the 1-d DOS 
p(e) = — e 2 ), using Eq.(6) leads to exactly the 

same function and prefactor as the exact solution (with 
U c = W/2 now being the H-I critical value). This sug- 
gests, that the H-I works particularly well in 1 dimension. 
Furthermore, exact results away from half- filling n = 1 
reveal, that there is a critical boundary separating states 
with doubly occupied sites from states with no such sites, 
which in the H-I framework considered here corresponds 
to the Lishitz lines Ul depicted in Fig.l. The exact 
boundary is given by the relation Ul = IFcos(7r(l — n)) 
[17]. Notice that, apart from the value of the Mott tran- 
sition point, the Lifshitz line obtained using the H-I ap- 
proach here is in good quantitative agreement with these 
results (see Fig.l). In general dimensions, the expected 
phase diagram is expected to be qualitatively similar [19] 



as in 1 dimension. We thus see, that the phase diagram 
derived by the Hubbard-I approach is in good qualitative 
agreement with exact results. 

Comparison with the Gutzwiller approach It is worth- 
while to compare our H-I results, for half filling, with 
those obtained using the Gutzwiller approximation to the 
model H c . We perform a standard calculation (see e.g. 
[12]), i.e assuming the reference state, that is Gutzwiller 
projected, to be the product of Fermi seas of the two 
spin species. Within the GA, the average energy of the 
Hamiltonian H c can be written as: 

(H c ) = q^oa + UD, (8) 

a 

where eo& is the energy of the Fermi sea of a given spin 
species. The band narrowing factor q a is easily found to 
be given by: 

_ (n a - D){1 + D - n a - n 9 ) + D(n s - D) 

1<y — T, x • l y J 

n a (l - n a ) 

Minimizing the average energy with respect to D, at 
half filling, we obtain a Gutzwiller-Brinkman-Rice tran- 
sition(GBR) at a finite U c = — 4 1 eo | , where eo is the 
summed ground state energy of the noninteracting Fermi 
liquids. The GA density of doubles on the metallic side 
is linearly dependent on the interaction U, reducing to 
zero at the transition point: 

D = ^(l-U/U c ), U<U C . (10) 

The GBR transition is of a quantitatively different 
character than the obtained H-I transition. Indeed, the 
distinct points are that the transition point is lattice de- 
pendent, while the metallic side has lattice independent 
double occupancies (which seem to be doubtful in the 
light of some exact results summarized above). As an 
aside, note however that, the choice of the reference state 
as noninteracting uncorrelated Fermi seas is of doubtful 
applicability, in the Gutzwiller analysis, of the model H c 
given that it certainly does not even correspond to the 
ground state of a correlated hopping Hamiltonian any- 
where, including in the limit [7 = 0. This is analogous 
to the Gutzwiller study of the bond-charge Hamiltonian 
Eq.(7) performed in [20], where the results for large bond 
charge interaction X cannot be considered reliable. In- 
deed, in particular at the symmetry point (which is of 
direct relevance to the model H c ), Xij = —Uj, Kollar 
and Vollhardt in [20] obtain a critical point U c — us- 
ing the GA. Thus, rather interestingly, we observe that 
the Gutzwiller approximation yields a qualitatively bet- 
ter picture of the phase transition itself in the model 
H c than in the bond-charge Hamiltonian. Finally, the 
GA shows a ferromagnetic instability in the bond-charge 
model, for certain lattices [20] before the Mott insulator 
transition. However, this is not the case for the GA in 
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the H c Hamiltonian. Indeed, as usual [6, 20] one can cal- 
culate the bulk magnetic susceptibility x, which is here 
given by: 



X 



1 - {U/Ucf 
2p(s F ) 



(1 - p(e F )U c ) 



(11) 



where p(ep) is the density of states at the Fermi sur- 
face. Only one factor depends on U and can diverge here, 
accompanying the metal insulator transition U —> U c . 
Thus, the GA, like the H-I calculation, does not describe 
a ferromagnetic metal before the Mott transition. How- 
ever, apart from this the phases and critical points are 
better described by the H-I approach. 

The Roth-corrected approach We have shown that the 
Hubbard-I approach to the model H c gives the correct 
scale of energy for the Mott transition. For general di- 
mensions, however, the critical point U c does not depend 
on magnetic polarization, since the gap is always given 
by U — W/2. Recall, that in obtaining the description of 
the Mott transition we neglected terms, in the equations 
of motion, related to magnetic order and other terms. 
All these terms can be treated in a unified manner by 
the procedure invoked by Roth [21], which removes am- 
biguities in the decoupling of Green's functions equations 
of motion method. Applying the Roth procedure at the 
Hubbard-I level of equations of motion for the consid- 
ered Hamiltonian H c , we obtain the following Green's 
functions: 



rg(w) 



1 n~ 



1 h* 



2-ku-E ' kK > 2ir u -E a 



(12) 



where, the quasi-particle energy factors are 



E a = (U - I-l) + ek{n os n Ss )/n s 
-€ k {{T D ) + (T x ))/n 9 + {{T UHB ) - (T LH B))/n 9 , (13) 

E a = -fi + ek{ho s hscj) /h s 
-e k ((T D ) + {T x ))/h ff + ({T UHB ) - (T LHB ))/h s . (14) 

The term (T x ) = (coa c 8& c 8cr c Ocr) describes the process 
of interchange of spins between neighbouring sites, while 
(Td) = (cba c S9 c o<j c 8cr) describes the process of doublon 
transfer. The averages in Eqs. (13, 14) are yet to be de- 
termined quantities, which is a standard feature of Roth's 
method. Notice, that disregarding the second lines in 
Eqs. (13, 14), and considering site occupations as uncor- 
rected, one recovers the Hubbard I results. 

A full analysis of the Roth solutions shall be considered 
elsewhere [22], while here we only consider implications 
for the Mott phase in which (T x ), (Td), (T UHB ), (Tlhb) 
vanish identically for the considered model H c . The clo- 
sure of the spectral gap to excitations indicates the Mott 
phase instability. Consider the case half-filled case with 
rif = nj, = 1/2. If there are no intersite correlations in 
the Mott states, U c = W/2 as indicated above. On the 
other hand for saturated ferromagnetic correlation, i.e. 



macroscopic separation of the system into two ferromag- 
netic domains, Eqs. (13, 14) show that there is no renor- 
malization of the band width in the Hubbard bands, and 
thus U c = W. As all Mott states have the same energy, 
this analysis indicates that the Mott phase is unstable 
already at U c = W (which is in good agreement with 
known results summarized earlier), and changes proba- 
bly into a ferromagnetic metallic state. 

Conclusions and Outlook In this work, we have an- 
alyzed the Hubbard I approximation, explicitly showing 
that its known drawbacks originate from the interband 
hopping T±i terms in the Hubbard model. We proposed 
to use the H-I approach in conjunction with perturba- 
tional expansion. The H-I approach, in combination with 
lowest order perturbational expansion, leads to a physi- 
cally appealing, picture of the Mott transition including 
the appearance of a extremely correlated Fermi liquid in 
its vicinity, which is complementary to the Gutzwiller- 
Brinkman-Rice picture. 

It is natural to wonder how well this approach fits as 
a description of Mott transitions and the surrounding 
ECFL in realistic strongly correlated electron models. In 
this regard, we emphasize that the double occupancy con- 
serving Hamiltonian H c analyzed here is equivalent to 
the Hubbard model with bond charge interactions at the 
symmetry point X = — t — which has been argued to be 
a quite realistic value (see e.g. [15, 17]). Our H-I or Roth 
improved H-I calculations compare quite favourably with 
known and expected results for the latter model, and are 
remarkably consistent with them in 1 dimension. Indeed, 
one would expect the picture presented here to hold in 
the neighbourhood of the symmetry point. On the other 
hand, for pure Hubbard like-systems, i.e. X close to 0, of 
course this strong coupling picture of the Mott transition 
can only be qualitatively correct (when antiferromagnetic 
order is suppressed) and that too only near the transition 
point. This gives rise to an open question concerning the 
possibility of obtaining an interpolating scheme between 
these two extreme cases. 

In this regard, we propose to use the following varia- 
tional ground state for systems exhibiting metal to Mott 
insulator transitions: 

|*> = exp[-afo]cxp[-77^]|* ) = exp[-af ]|GWF), 

(15) 

where a,n arc variational parameters and l^o) is an 
appropriate reference state. This state is an extension 
of the standard Gutzwiller Wave Function (GWF) and 
should provide an improved description of the metal- 
lic side of the transition. Note that the two exponen- 
tial terms commute and together form the exponential 
of —aH c (U = rj/a), which may be viewed as a partial 
projection on to the ground state of H C (U = tj/a). This 
makes a connection with the ECFL properties described 
in this paper and allows for dressing of the GWF with 
precursors of the Hubbard bands. It is worth mentioning 
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here that a similar function (containing hopping only of 
the lower Hubbard band in To) has already been used 
in [23] and has provided excellent results in comparison 
with exact results for 2 electrons on the Hubbard square. 
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